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.Abstract 

■ Magnetic excitations and spin correlations near the interface of two spin-^ Heisenberg antiferromagnets are considered 
|using the spin-wave approximation. When the interaction between boundary spins differs essentiahy from exchange 
■constants inside the antiferromagnets, quasi-two-dimensional spin waves appear in the near-boundary region. They eject 
[bulk magnons from this region, thereby dividing the antiferromagnets into areas with different magnetic excitations. 
■The decreased dimensionality of the near-boundary modes leads to amplified nearest-neighbor spin correlations in the 
[interface area. 

[Keywords: Heterostructure, semi-infinite Heisenberg antiferromagnet, magnetic excitations, spin correlations 
■PACS: 75.10.Jm, 75.30.Ds, 75.70.Cn 
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1. Introduction 

In recent years, an active interest has been taken in het- 
erostructures fabricated out of magnetic crystals. Looking 
for new effects and their possible applications, a wide vari- 
ety of systems has been investigated both experimentally 
and theoretically (see, e.g., [l|-l5|)- One of the systems 
belonging to this group is the interface of a Heisenberg 
antiferromagnet with vacuum. Spin excitations and corre- 
lations near this interface were investigated in Refs. 
In particular, it was shown that the antiferromagnet is di- 
vided into two regions with different spin excitations 
ll| . In the first region, which embraces two near-boundary 
atomic layers, the excitations have dimensionality one less 
than the crystal dimensionality. The rest of the crystal is 
the second region, which elementary excitations are stand- 
ing spin waves. The near-boundary excitations eject these 
latter spin waves from the near-boundary region, which 
leads to the mentioned separation of the antiferromagnet 
into the two regions. Both physically and mathematically 
this situation resembles the problem of the local impurity 
considered by I.M. Lifshits [lH, in which localized elec- 
tronic states eject bulk states from the impurity region. 
The decreased dimensionality of the near-boundary modes 
leads to increased spin correlations in the interface region 
and their damped oscillations away from the region [6r[8|. 

In this Letter the heterostructure of two semi-infinite 
spin-i Heisenberg antiferromagnets with nearest-neighbor 
exchange constants Ji and J2 = O.8J1 is considered. Such 
heterostructures were already fabricated 13|, however, 
to our knowledge magnetic excitations near the interface 
were not experimentally investigated yet. Both antiferro- 
magnets are supposed to have identical simple cubic lat- 
tices with the interface coinciding with a crystallographic 



plane. As in the case of the antiferromagnet- vacuum inter- 
face, for the exchange constant between sites of the two an- 
tiferromagnets J3 <C J2 modes of decreased dimensionality 
exist in the near-boundary region, expelling bulk modes 
from this area. These near-boundary modes are seen as 
sharp peaks in the spectral function characterized by the 
layer index and the wave vector k describing movement 
parallel to the interface. With J3 approaching J2 the near- 
boundary mode disappears at first in the antiferromagnet 
with J — J2 and for small k. With further growth of J3 
the mode fades away also in the second antiferromagnet. A 
new type of near-boundary excitations appears at J3 « Ji 
in the antiferromagnet with J = J2, and at a somewhat 
larger J3 in the other antiferromagnet. Their frequencies 
surpass frequencies of bulk spin waves. These interface ex- 
citations are mainly connected with pairs of spins on each 
side of the boundary. In the near-boundary region, the 
nearest-neighbor spin correlations are amplified in com- 
parison with the bulk value, when J3 <^ J2 and J3 > Ji. 
The deviations nearly disappear when J3 approaches J2. 
At J3 ~ (Ji -I- J2)/2 the deviations reverse their signs. 



2. Main formulas 

The axes are chosen in such a way that the antiferro- 
magnet with the exchange constant Ji is located in the 
half-space Ix ^ 0, while the other antiferromagnet with 
the exchange constant J2 is in the half-space Ix ^ — 1. 
Additionally there is an exchange interaction J3 between 
boundary spins of the antiferromagnets. Here sites of a 
three-dimensional simple cubic lattice are labeled by the 
three coordinates lx,ly,lz, and the lattice spacing is set as 
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the unit of length. The system is described by the Hamil- 
tonian 



H 



■h 



la /,>0 



J: 



la 



J3 X ^10^1,-1 , 



(1) 



where 1 = {ly, Iz), a — (±1, 0), (0, ±1) are four unit vectors, 
which connect nearest neighbor sites in the YZ plane, and 
Sl is the spin-i operator. 

Since for low temperatures the system has the long- 
range antiferromagnetic order, its low-lying elementary ex- 
citations can be described in the spin wave approximation, 

5£ = e'"^Q-6i,6L), S^^Pib^ + P^bl (2) 

where the spin- wave operators 6l and 6^ satisfy the Boson 
commutation relations, and 



n = (7r,^,7r), P± = - (1 ± e^n^) , L = {l^jyj,). 

If we substitute Eq. ^ into Eq. ([T]), drop constant terms 
and terms containing more than two spin- wave operators, 
we obtain 



H — Hi + H2 + -ffs, 
Hi = 3Ji J2 (^-Wo)blj 



t ;,t 

-k,L 



Ji 
2 



H2 = 3J2 J2 



kJ^sS-l 



(3) 



J2 E lkHb^'J-Kl.+bl^b\i^ 



k/^^-l 



k/x^-l 



^3 = YE(^ko^ko + &Li6k.-i 



Here we took into account the translational invariance of 
Hamiltonian (jlj in the YZ plane and used the Fourier 
transformation 



where k is a two-dimensional (2D) wave vector and N is 
the number of sites in the periodic YZ plane. In Eq. ([3]), 
^k^ = l[cos(fcy) + cos(fc^)]. 

To investigate the spectrum of elementary excitations, 
we introduce the two-component operator 



kU 



bki, 
5^ 



and define the matrix retarded Green's function 



(4) 



where Ek;,(t) = e'"'B^i^e 



Ht 



with H determined by 



Eq. ®. 

To calculate Green's function (|4]), we use the equation 
of motion 



ri 



3f3 + 27(') 



r2TiD{kt, I, + 1, Q + rshDO^t, I, - 1, Q 



J3-J1 



J3-J2 



J3 



2 ■ - 2 



where 

n 

Tl 



J20{-l^-l) + Jie{l,), 

i [J2e{-i, - 2) + j^eii,)] . 



(5) 



[J2^(-?x - 1) + JlOil:, - 1)] . 
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-1 
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Equation ([5|) contains local perturbations proportional to 
Si^o and By analogy with Lifshits' solution 

Green's function ^ can be expressed through the func- 
tion D'-^^ (ktlxl'-i.) of the equation without the perturbation 
terms. After the Fourier transformation and some mathe- 
matical manipulations this expression reads 

D(ku:lX) ^ D^'^HkulX) 
J3-J1 



2 

J3-J2 



D^"^kujl^O)D{kujOQ 
D("\kLul,,-l)D{kLu,^l,Q 



' 2 



-D^°\kujl^,-l)T2D(kujOQ, 
where fo is the 2x2 identity matrix and 



(6) 



T2 = 



1 

1 



^i?L)(")(kwOO)f2 



-1 



i? = ^(")(kw,-l,-l)f2 



To - 

i)(°)(kwOO). 



"^3 - -/l f^(0) 



D^">{ku;00) 



To calculate Green's function l)'^''^ (kw/^^i) it is neces- 
sary to diagonalize Hamiltonian ([3]) without the terms pro- 
portional to J3, Si^o, and This can be fulfilled us- 
ing the Bogoliubov-Tyablikov transformation 15|. The 
obtained Green's function D'^^'> (kujlxl'^) reads 11| 



0iQ0iQ Jdk, sin[fc,(?, + l)] 


X sin[fc,(C-|-l)]rkfe,(Ji), 

D^"^\kwlJ'J - - l)OH'^ - 1) y dk, sm{kJx) 





(7) 



Qkk^{J) 



XJ) 



- Ekk, (J) +11] UJ + Ekk^ (J) -I- iri 

AvkAJ)BkkAJ) Bl,JJ) 

BLM) A^kAJ)BwkAJ) 

AkkAJ)BkkAJ) Al^U) 



AkkAJ) = 



BkkAJ) 



Ekk^ (J) 



where 7^'^^ ~ -[cos(fc:r) + cos(fcj,) -I- cos(A:z)] and rj = +0. 
^ 3 
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Figure 1: The spectral function lmDii{]<.Lulxlx) for (a) J3 = 0.2, 
(b) J3 = 0.7 and (c) J3 = 1.1. k = (0.67r,0) and J2 = 0.8. Red 
solid and blue dash-dotted lines correspond to Ix = and Ix = 
— 1, respectively. In part (a), green dashed and purple doted lines 
show spectral functions for Ix = ^ and 2. In part (c), pink dotted 
and cyan dashed lines demonstrate spectral functions of unbounded 
antiferromagnets with J = Ji and J = J2, respectively. 



3. Spin excitations 

Below the exchange constant Ji is set as the unit of en- 
ergy, and the case J2 = 0.8 is considered. For J3 = 
the problem reduces to two noninteracting semi-infinite 
antiferromagnets , which spin excitations were considered 
in Refs. 0, 



Hi 



As mentioned in the Introduction, the 
spectrum of each antiferromagnet consists of a quasi-2D 
mode of spin waves, which are mainly located in two near- 
boundary layers, and bulk modes of standing spin waves 
in the rest of the crystal. This picture retains for J3 ^ J2. 
This case is illustrated in Fig. [Ua). In the boundary lay- 
ers Ix = and = — 1 of both semi- infinite antiferro- 
magnets the spectral function lmDii{kujlxlx) has sharp 
peaks, which rapidly lose intensity with distance from the 
interface, as seen from the spectral functions for — 1 
and 2 [the functions for — —2 and —3, which are not 
shown in Fig. [Ija), look analogously with the respective 
frequency shift]. These peaks correspond to the near- 
boundary modes, which propagate along the interface with 
the dispersion similar to that of the 2D magnons, as seen in 
Fig. [21 In the boundary layers the continuum of standing 
spin waves is seen as weak shoulders on the high-frequency 
sides of the peaks. However, already in layers Ix = 2 and 
Ix = —3 the peaks become weak and spectral functions be- 
come close to those in unbounded antiferromagnets shown 
in Fig. [TJc) [maxima in these functions are connected with 
power divergences in luiDnikojlxlx) at the edges of the 
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Figure 2: The dispersion of the near-boundary modes along the sym- 
metry Knes of the 2D Brillouin zone in the antiferromagnet with 
J = Jl (the black solid line) and in the antiferromagnet with J = J2 
(the blue dashed line). The red dash-dotted line demonstrates the 
dispersion of the 2D spin waves for J = J\. J3 = 0.2 and J2 = 0.8. 
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Figure 3: The spectral function Im Dn (kwiaj/i,) for k = (O.Svr, 0) 
(a,c) and for k = (0.67r, 0) (b,d). J3 = 1.51, J2 = 0.8. Red solid 
lines correspond to Ij: = and Ix = —1, blue dashed lines to Z^, = 1 
and Ij: = —2 on the lower and upper panels, respectively. Pink dotted 
lines are spectral functions of an unbounded antiferromagnets with 
J = Jl and J = J2 . 



spectra i?kfe^(>/); a finite artificial broadening transforms 
the divergences into maxima]. Tlius, for J3 ^ J2 the crys- 
tals are split into two regions with different spin excitations 
- four layers closest to the interface are the domain of ex- 
istence of the near-boundary modes with the suppressed 
intensity of bulk modes, while these latter modes dominate 
in the rest of the crystals. 

With increasing J3 the maximum of the near-boundary 
mode in the spectral function weakens and disappears 
starting from the antiferromagnet with the smaller ex- 
change constant and small wave vectors k. The near- 
boundary mode fades away at J3 ~ 0.5 in the antifer- 
romagnet with J = J2 and at J3 « 0.65 in the other 
antiferromagnet for all wave vectors. For J3 « J2 even in 
the near-boundary layers the spectrum consists of a contin- 
uum of bulk states [see Fig.[TIb)]. New peaks start to form 
near the upper edges of the continuum spectra at J3 « Ji 
in the antiferromagnet with J — J2, and at J3 « 1.2 in 
the second antiferromagnet [see Fig.[TJc)]. As follows from 
Fig. 3, these peaks are seen only in the interface region 
and, therefore, they correspond to near-boundary modes. 
For moderate J3 only one of the two peaks and its replica 
in the neighbor antiferromagnet may be observed for some 
momenta k [see Figs.|3Ua) and (c)]. For other k both near- 
boundary peaks and their replicas are seen [Figs.[2Ib) and 
(d)]. In the antiferromagnet with J ~ Ji the dispersion 
of the mode resembles in shape the dispersion of the 2D 
spin waves, however, with a finite offset of frequencies by 
Wo « 3Ji, 



Wk = c^o + 2J'yi- (7^'') . (8) 

For parameters of Fig. 3 J' w 0.13, much smaller than Ji. 
In this respect this mode differs from the near-boundary 
modes existing at J3 <^ J2, which dispersions correspond 



to effective exchange constants larger than Ji (see Fig. [5]). 
Apparently in the limit of very large J3 the near-boundary 
modes in Fig.[2]transform to triplet excitations of the pairs 
of boundary spins. 

4. Spin correlations 

Spin correlations between nearest neighbors can be ex- 
pressed in terms of correlations of spin- wave operators us- 
ing Eq. ^ and the translational invariance of Hamiltonian 
^ in the YZ plane, 

(SlSlO = ^^|2cos[k(l-l')](6ki.6-k,ii> 

k 

+ {blM^) + {bl,b^iO]-\- (9) 
Bearing in mind the property of Green's function ^ 

the spin-wave correlations in Eq. (|9]) can be expressed as 

{B^i.B^i-} - , (10) 

where f3 — T^^ is the inverse temperature. 

Let us denote the spin correlations on sites located par- 
allel and perpendicular to the interface as 

Cl{Ix) = (Sl+a,i^Si/_^) , 

Ct (^Ix + 2) ^ (Si,;,+iSi/_^) . 
For T = with the use of Eq. (fTU|) these spin correlations 



4 




„. ^ mi I ■ , ■ r ,1 , ■ , , ■ I to ACt(— 0.5), the left y axis to all other quantities, 

figure 4: ihe deviations oi the nearest-neighbor spin correlations 

from the bulk value as functions of Ix for sites located parallel ( ACl , 

a) and perpendicular (ACy, b) to the interface. J2 = 0.8, values of 

J3 are indicated in the figure. 
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Let us define deviations of spin correlations parallel and 
perpendicular to the interface from the bulk spin correla- 
tion Cb = -0.3005 [i, [13 as 



Cl,t — Cb 

Cb ■ 



The calculated values of these deviations are shown in 
Figs, m and [S] for T = 0. As seen from these figures, main 
deviations of the spin correlations from the bulk value fall 
on the boundary and the second to the boundary layers, 
i.e. on the existence domain of the near-boundary modes. 
The largest in absolute value deviations are obtained in 
and between the boundary layers for J3 ^ J2 and for 



J3 > Ji, when the near-boundary modes are well resolved 
in the spectral function. As indicated in Refs. [ol- 11 . 
in the case of a semi-infinite antiferromagnet the ampli- 
fied spin correlations in the boundary layer are connected 
with the quasi-two-dimensionality of the near-boundary 
mode, which dominates in the spectral function of this 
layer. It behaves like a 2D antiferromagnet, in which the 
nearest-neighbor spin correlation is larger in modulus than 
in a three-dimensional antiferromagnet (0.3346 vs. 0.3005 
16 . l3|)- This explains the strengthened spin correlations 



in this layer and between it and the second to the bound- 
ary layer. The correlations in the second and between the 
second and third layers are smaller in modulus than Cb due 
to the destructive contribution of the boundary and bulk 
modes. For the case J3 ^ J2 in both antiferromagnets 
the distribution of spin correlations is similar to that men- 
tioned above and apparently has the same explanation. 
Peaks of the near-boundary mode disappear at J3 w 0.5 
in the antiferromagnet with J = J2 and at J3 ~ 0.65 
in the other antiferromagnet, which allows bulk modes to 
penetrate in the near-boundary region. As a consequence 
spin correlations in the second and between the second and 
third layers of an antiferromagnet become approximately 
equal to Cb- The new type of interface excitations appears 
at J3 « Ji in the antiferromagnet with J — J2 and at 
J3 « 1.2 in the other antiferromagnet. In this mode, spin 
correlations between sites on each side of the interface and 
in the second and between the second and the third lay- 
ers are amplified, while correlations on the boundary and 
between the boundary and the second layers are weakened 
in both antiferromagnets. As follows from Figs. |4] and [5l 
the deviations of the spin correlations from the bulk value 
change sign at J3 « ( Ji + J2)/2. 
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5. Conclusion 

In this Letter, magnetic excitations and nearest- 
neighbor spin correlations near the interface of two semi- 
infinite spin-i Heisenberg antiferromagnets were inves- 
tigated using the spin-wave approximation. It was as- 
sumed that the antiferromagnets have identical simple cu- 
bic lattices, the nearest-neighbor exchange constants Ji 
and J2 = 0.8 Ji in the bulk of the crystals and J3 between 
the interface sites. It was shown that for J3 ^ J2 a near- 
boundary region arises, which embraces four layers clos- 
est to the interface. In this region, magnetic excitations 
are quasi- two-dimensional spin waves, which expel bulk 
modes from their domain. As J3 approachcis J2 peaks of 
the near-boundary modes lose intensity and finally disap- 
pear in the continuum of the bulk modes, which penetrate 
in the near-boundary region. This process starts in the 
antiferromagnet with the smaller exchange constant and 
at small wave vectors. At J3 ~ Ji in the antiferromagnet 
with J = J2 and then, at J;^ w 1.2Ji in the other anti- 
ferromagnet a new type of near-boundary modes arises. 
In these modes pairs of spins on each side of the interface 
are mainly involved. The separation of the heterostructure 
into the regions with diflFerent spin excitations leads to a 
peculiar distribution of nearest-neighbor spin correlations 
near the interface. In the case J3 <C J2 the decreased 
dimensionality of the near-boundary modes leads to the 
correlations in the boundary layers, which are increased 
in modulus in comparison with the bulk correlation. At 
the same time the destructive contribution of the near- 
boundary and bulk modes decreases correlations between 
some sites. For J3 « (Ji + J2)/2 the correlations are close 
to the bulk one. In the case J3 > Ji the correlations are 
amplified between the interface spins, and they are weak- 
ened in the boundary layers of both antiferromagnets. The 
deviations of the correlations from the bulk value reverse 
sign at J3 w (Ji -I- Ji)/2. 
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